Varying- Alpha Cosmologies with Potentials 
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We generalize the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) model for the variation of the 
fine structure 'constant', a, to include an exponential or inverse power-law self-potential for the 
scalar field tp which drives the time variation of a, and consider the dynamics of tp in such models. 
We find solutions for the evolution of <p or a in matter-, radiation- and dark-energy-dominated 
cosmic eras. In general, the evolution of tp is well determined solely by either the self-potential or 
the coupling to matter, depending on the model parameters. The results are general and applicable 
to other models where the evolution of a scalar field is governed by a matter coupling and a self- 
potential. We find that the existing astronomical data stringently constrains the possible evolution 
of a between redshifts z ~ 1 — 3.5 and the present, and this leads to very strong limit on the allowed 
deviation of the potential from that of a pure cosmological constant. 

PACS numbers: 98.80.Es, 98.80.Bp, 98.80.Cq 



I. INTRODUCTION 

For the first time there is a body of detailed astronomi- 
cal evidence consistent with the time variation of a tradi- 
tional constant of Nature. The observational programme 
of Webb et al. [l|, Q has completed detailed analyses of 
128 Keck-HIRES quasar absorption line systems at red- 
shifts 0.5 < z < 3 using the many-multiplet method to 
compare separations between line separations affected by 
special-relativistic effects and found evidence consistent 
with the fine structure constant at redshift z, a(z), hav- 
ing been smaller in the past, at redshifts z ~ 1 — 3.5. The 
shift in the value of a between its value a(z) at redshift 
z and its present-day value cv(0), for all the data is given 
provisionally by 

Aa/a ee [a{z) - a(0)]/a(0) = (-0.57 ± 0.10) x 10~ 5 . 

Subsequent reduction of another data set of 23 VLT- 
UVES quasar absorption systems 0.4 < z < 2.3 by 
Chand et al, [1, HJ using a partial version of the many- 
multiplet method at first produced a result consistent 
with no variation in a, with an unusually small uncer- 
tainty, Aa/a = (-0.06 ±0.06) x 10" 5 . However, the data 
reduction did not allow Aev to be a free parameter in the 
data fitting, and a reanalysis of the same data set by 
Murphy et al. @] using the full many-multiplet method 
increases the uncertainties sixfold, and leads to a revised 
bound of 

Aa/a = (-0.64 ± 0.36) x 10~ 5 . 

Any present-day variation of a can also be constrained by 
direct laboratory comparisons of clocks based on different 
atomic frequency standards over a period of months or 
years. Until recently, the most stringent atomic clock 



constraints on any current temporal variation of a were 

a/a = (-3.3 ± 3.0) x 10~ 16 yr' 1 , 

which arose by combining measurements of the frequen- 
cies of Sr 0, Hg 0], Yb 01, and H [| relative to Cae- 
sium; Cingoz et al. [l(| have also recently reported a 
less stringent limit of a/a — — (2.7± 2.6) x 10~ 15 r/r~ 1 . 
If the systematic errors can be fully understood, an ul- 
timate sensitivity of 10~ 18 ?/r~ 1 may be possible with 
this method fill ]. If a linear variation in a is assumed 
then the Murphy et al. quasar measurements equate to 
a/a = (6.4 ± 1.4) x lO" 16 ?^- 1 [J, 0. If the variation 
is due to a light scalar field described by a theory like 
that of Bekenstein and Sandvik, Barrow and Magueijo 
(BSBM) [13, [lj], then the rate of change in the con- 
stants is exponentially damped during the recent dark- 
energy-dominated era of accelerated expansion, and one 
typically predicts a present-day value of 



a/a = 1.1 ±0.3 x 10" 



yr 
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by direct extrapolation from the Murphy et al. data 
d This is not ruled out by the atomic clock con- 
straints mentioned above. For comparison, the Oklo nat- 
ural reactor constraints, which are based on the need for 
the Sm 149 + n — > Sm 147 + 7 neutron capture resonance 
at 97.3 MeV to have been present 1.8 — 2 Gy r ago at 
z = 0.15, as first pointed out by Shlyakhter [14]], are cur- 
rently [l| Aa/a = (-0.8±1.0)xl0 l8 or(8.8±0.7)xl0- 8 
(because of the double-valued character of the neutron 
capture cross-section with reactor temperature) and [l6| 
Aa/a > 4.5 x 10~ 8 (6a), when the non-thermal neu- 
tron spectrum is taken into account. However, there 
remain significant environmental uncertainties regarding 
the reactor's early history and the relationship between 
changes in the resonance energy level and those in the 
values of any underlying constants. For reviews of the 
wider issue of varying constants in addition to a, see the 
reviews in refs. jl7J , and for some implications of the uni- 
fication of fundamental forces see refs. Hal. 
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Recently, Rosenband et al. [19| measured the ratio of 
aluminium and mercury single-ion optical clock frequen- 
cies, /Ai+//H g +, at intervals over a period of about a 
year. From these measurements, the linear rate of change 
in this ratio was found to be (—5.3 ± 7.9) x 10 -17 yr _1 
(but see ref. for some refinements). These measure- 
ments provides the strongest limit yet on any temporal 
drift in the value of a: 

a/a = (-1.6 ± 2.3) x KT 17 yr' 1 . 

This limit is strong enough to exclude theoretical expla- 
nations of the change in a reported by Webb et al. [l], 0] 
based on the slow variation of an effectively massless 
scalar field [13], even allowing for the damping by cos- 
mological acceleration, unless there is a significant effect 
that slows the locally observed effects of changing a on 
cosmological scales (for a detailed analysis of glob al-local 
coupling of variations in constants (see Refs. [2l|). 

Theories in which a varies will in general lead to vio- 
lations of the weak equivalence principle (WEP). This is 
because the a variation is carried by a scalar field, (p, and 
this couples differently to different nuclei because they 
contain different numbers of electrically charged particles 
(protons). The theory discussed here has the interesting 
consequence of leading to a relative acceleration of order 
10~ 13 [22] if the free coupling parameter is fixed to the 
value given in Eq. (O using a best fit of the theories cos- 
mological model to the quasar observations of refs. [l], 0] • 
Other predictions of WEP violations have also been made 
in refs. [H, [24], [25[ . The observational upper bound on 
this parameter from direct experiment is just an order of 
magnitude larger, at 10 -12 , and limits from the motion of 
the Moon are of similar order, [26[, but space-based tests 
planned for the STEP mission are expected to achieve a 
sensitivity of order 10 -18 and will provide a completely 
independent check on theories of time-varying e and a 
[221123. 

In view of this tension between direct local measure- 
ments and astronomical measurements of the fine struc- 
ture 'constant' it is important to explore the widest 
possible range of self-consistent theoretical models for 
the time-evolution of a so as to understand the pos- 
sible evolutions of Aa/a over the range < z < 6 
that spans the astronomical, geochemical, and labora- 
tory measurements. In the remainder of this paper we 
will present cosmological extensions to the simple BSBM 
scalar field models for varying a that include a non- 
zero self-interaction potential, V(ip) for the scalar field, 
<p, carrying the spacetime evolution of a. We will con- 
sider two representative theories, where V has exponen- 
tial and power-law variation, respectively, and determine 
the solutions for the cosmological evolution and the time- 
variation of a during the radiation, dust, and dark-energy 
dominated eras of the universe. 

The organization of this paper is as follows: in § Inland 
§ 11111 we present the theory of varying a based on the cou- 
pling of a scalar field to the electromagnetically charged 
matter and list the relativistic equations for the investi- 



gations of this theory. In § IIVI and §[V]we review, respec- 
tively, the cosmological evolutions of the scalar field tp for 
the model with no scalar field self-interaction potential 
(just coupling with matter) and for quintessence models 
with exponential and inverse power-law self-potentials. 
The § IVII is devoted to an investigation of how the scalar 
field (p evolves if both the matter coupling term and the 
bare self-potential are non-zero, which is supplemented 
by the numerical examples shown in § I VIII Finally, con- 
clusions are drawn in § I Villi 

II. BSBM SCALAR-FIELD THEORIES FOR 
VARYING a 

There are a number of possible theories allowing for 
the variation of the fine structure constant, a. In the 
simplest cases we take c and H to be constants and at- 
tribute variations in a to changes in e or the permittiv- 
ity of free space (see [29j] for a discussion of the meaning 
of this choice). This is done by letting e take on the 
value of a real scalar field which varies in space and time. 
Thus eo — ► e = e^eix^), where e is a dimensionless scalar 
field and eo is a constant denoting the present value of 
e. This operation implies that some well established as- 
sumptions, like charge conservation, must give way [30I ]. 
Nevertheless, the principles of local gauge invariance and 
causality are maintained, as is the scale invariance of the 
e field (under a suitable choice of dynamics) and there is 
no conflict with local Lorentz invariance or covariance. 

The dynamics are then constructed as follows. Since 
e is the electromagnetic coupling, the e field couples 
to the gauge field as eA^ in the Lagrangian and the 
gauge transformation which leaves the action invariant is 
eA^ — > eA^ + rather than the usual A^ — ► A^ + x,n- 
The gauge-invariant electromagnetic field tensor is there- 
fore 

F„u = \ {{eA v ) tll - (eA„) iU ) , (1) 

which reduces to the usual form when e is constant. The 
electromagnetic part of the action is still 

S em = ~J d^x^—gF^F^. (2) 

and the dynamics of the e field are controlled by the 
kinetic term 

S. = ~^JdW=g^, (3) 

as in dilaton theories. Here, I is the characteristic length 
scale of the theory, introduced for dimensional reasons. 
This constant length scale gives the scale down to which 
the electric field around a point charge is accurately 
Coulombic. The corresponding energy scale, hc/l, has 
to lie between a few tens of MeV and the Planck scale, 
~ 10 19 GeV to avoid conflict with experiment. This 
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generalisation of the scalar theory proposed by Beken- 
stein | 1 2 fl was made by Sandvik, Magueijo and Barrow 
|SE [3JI3I [ai| and will be referred to as the BSBM the- 
ory. It includes the gravitational effects of ip and gives 
the field equations: 

= SirG (T M T tCr + T* v + T^e^) . (4) 

The stress tensor of the <p field is derived from the la- 
grangian C v — — ^d^ipd^tp and the ip field obeys the 
equation of motion 

Dip = -e- 2 ^C cm (5) 
u 

where we have defined the coupling constant uj — (c)/l 2 . 
This constant is of order ~ 1 if, as in [l3j], the energy 
scale is similar to the Planck scale. It is clear that £ em 
vanishes for a sea of pure radiation since then £ om = 
(E 2 — B 2 )/2 — 0. We therefore expect the variation in a 
to be driven by electrostatic and magnetostatic energy- 
components rather than electromagnetic radiation and 
with h = c = 1, the fine-structure 'constant' is given by 

a/ao = e 2 /eQ = exp(2ip). 

The considerations raised by Duff [35[ do not impact 
upon well-defined varying 'constant' theories like this, 
even if they appear dimensionful. The presence of a new 
field, like ip, always requires a second-order energy con- 
servation equation, like Eq. ([5]) and the integration of this 
equation always leads to a new integration constant, tpo, 
with the same dimensions as tp and so the evolution of 
the dimcnsionless quantity <p/<po involves no ambiguities 
under redefinitions of units. 

In order to make quantitative predictions we need 
to know how much of the non-relativistic matter con- 
tributes to the RHS of Eq. ([5]) . This is parametrised by 
£ = C em /p, where p is the energy density, and for bary- 
onic matter C em = E 2 /2. For protons and neutrons £ p 
and can be estimated from the electromagnetic cor- 
rections to the nucleon mass, 0.63 MeV and —0.13 MeV, 
respectively [24|. This correction contains the E 2 /2 con- 
tribution (always positive), but also terms of the form 
j^,a^ (where j M is the quarks' current) and so cannot 
be used directly. Hence, we take a representative value 
Cp ~ Cn ~ 10~ 4 . Furthermore, the cosmological value 
of £ (denoted £ m ) has to be weighted by the fraction of 
matter that is non-baryonic. Hence, ( m depends strongly 
on the nature of the dark matter and can take both posi- 
tive and negative values depending on which of Coulomb- 
energy or magnetostatic energy dominates the dark mat- 
ter of the Universe. It could be that (cdm « — 1 (super- 
conducting cosmic strings, for which C em ~ — i? 2 /2), or 
Ccdm *C 1 (neutrinos). BBN predicts an approximate 
value for the baryon density of fls w 0.03 (where fie 
is the density of matter in units of the critical density 
3H 2 /8irG) with a Hubble parameter of H = 60 Kms" 1 
Mpc^ 1 , implying flcDM ~ 0.3. Thus, depending on the 
nature of the dark matter, £ m can be virtually anything 



between —1 and +1. The uncertainties in the underlying 
quark physics and especially the constituents of the dark 
matter make it difficult to impose more certain bounds 

on Cm- 

There are a number of conclusions that can be drawn 
from the study of the simple BSBM models with ( m < 0. 
These models gave a good fit to the varying a implied 
by the quasar data of refs. [l|, |2j. There is just a single 
parameter to fit and this is given by the choice fl3j 

- ^ = (2±1) x 10~ 4 (6) 

The simple solutions of the BSBM theory predict a slow 
(logarithmic) time increase of a during the dust era of 
k = Friedmann universes. The cosmological constant 
turns off the time- variation of a at the redshift when the 
universe begins to accelerate (z ~ 0.7) and so there is 
no conflict between the a variation seen in quasars at 
z 1 — 3.5 and the limits on possible variation of a 
deduced from the operation of the Oklo natural reactor 
[3 EH (even assuming that the cosmological variation 
applies unchanged to the terrestrial environment). The 
reactor operated 1.8 billion years ago at a redshift of only 
z ~ 0.1 when no significant variations were occurring in 
a. The slow logarithmic increase in a also means that 
we would not expect to have seen any effect yet in the 
anisotropy of the microwave backgrounds [ly, H3, l38j : 
the value of a at the last scattering redshift, z = 1000, 
is only 0.005% lower than its value today. Similarly, the 
essentially constant evolution of a predicted during the 
radiation era leads us to expect no measurable effects on 
the products of Big Bang Nucleosynthesis (BBN) [3!| be- 
cause a was only 0.007% smaller at BBN than it is today. 
This does not rule out the possibility that unification ef- 
fects in a more general theory might require variations in 
weak and strong couplings, or their contributions to the 
neutron-proton mass difference, which might produce ob- 
servable differences in light-element nucleosynthesis, and 
new constraints on varying a, at z ~ 10 9 — 10 10 . By 
contrast, varying-alpha cosmologies with £ > lead to 
bad consequences unless the scalar field driving the al- 
pha variations is a 'ghost' field, with negatively coupled 
kinetic energy, in which case there can be interesting cos- 
mological consequences, [4(|. The fine structure 'con- 
stant' falls rapidly at late times and the variation is such 
that it comes to dominate the Friedmann equation for the 
cosmological dynamics. We regard this as a signal that 
such models are astrophysically ruled out and perhaps 
are also mathematically badly behaved. 

The earlier analyses of the cosmological solutions of the 
BSBM equations considered the situation in which the 
scalar field driving variations in a has no self-interaction 
potential, V(tp) = 0. In this paper, we are going to ex- 
plore some of the consequences for the time- variation of 
a = exp(2ip) that arise when we introduce a non-zero 
potential for the scalar field driving the variations in a. 
We will consider the representative cases of the exponen- 
tial potential V(tp) = Vq exp(— \<p) and inverse power- 



4 



law potential V[tp) — V^tp^ 7 , and and classify the new 
behaviours that arise for A, 7 ^ 0. We note that the cases 
of A, 7 = correspond to V = Vq, which is equivalent to 
the presence of a cosmological constant. The solutions 
for such scenarios were found in our earlier studies and ip 
relaxes quickly to a constant asymptotic value once the 
expansion starts to accelerate. 



III. COSMOLOGICAL EQUATIONS 

In the BSBM theory, the total action of the Universe 
is given by 



S = 



-<?(£ 



grav 



e 



-2.p 



The universe is described by a homogeneous and isotropic 
Friedmann metric with expansion scale factor a(t). The 
Friedmann equation is given by 

H 2 = \ [ Pm (l + \C\e- 2 n + Pr e- 2 * + Pv + pa] ; (7) 

where we assume that the universe is spatially flat; the 
quantities p m , p r , p v , p\ are the energy densities in non- 
relativistic matter, relativistic matter, scalar field and 
cosmological constant (so p\ is a constant), respectively. 
We will first consider the case where the scalar field tp 
has no potential term, and then consider the cases with 
exponential and power-law potentials. 

The conservation equations for matter and radiation 
are given as 



p m + 3Hp m = 0, 
p r + AHp r = 2p r tp, 

and the scalar field equation of motion is 

■■ , orr . , dV(p) 2|C| 

tp + 3Htp+ — = p m e v . 

otp UJ 



(8) 
(9) 



(10) 



Taking the time derivative of Eq. (7), and using 
Eqs. (8)-(10), we get 



H = — 
2 



Pm 



1 + \(\e- 2 *) + ^p r e~ 2 * -! ^ 



.(11) 



In a universe filled with non-relativistic matter and the 
scalar field, p\ = Pr = this equation reduces to 



V 2 ~ Vdp) 



IV. THE CASE OF ( < 0: V CONSTANT 

This was the situation analysed in the original presen- 
tation of the BSBM theory in Refs. [U, I32TI331. 134] ] and 
extended to include higher-order corrections in Ref. [43| , 



small perturbations [4l| , and a linearised potential in ref. 
[42l ]. The structure of the cosmological solutions has an 
expected feature. The cosmological dynamics of the scale 
factor a(t), controlled by the Friedmann equation, is not 
influenced to leading order by the small variations in tp. 
However, the cosmological variation of a(t) has a sig- 
nificant effect on the dynamics of ip, and hence upon 
the evolution of the fine structure 'constant' a(t). The 
key results for a cosmological model that evolves through 
a radiation-CDM-vacuum energy-dominated sequence of 
three phases are as follows: 

During the radiation era in which a — t x l 2 , there is an 
exact solution of 

(tpa 3 ) = Nexp{-2tp) 
where N > is a constant defined by 

N = p m a A 

ijj 



given by 



<p = \ log(87V) + i log(t) 



For physically realistic choices of the parameters the log- 
arithmic term is never significant during the radiation 
era of our universe and tp is constant then, which would 
be expected since E 2 = B 2 for the radiation equilibrium 
which means that ( ~ (E 2 - B 2 )/(E 2 + B 2 ) is effectively 
zero and ip constant. 

During the matter-dominated era a = i 2 / 3 , and there 
is a late-time asymptotic series solution of the form 
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1 1 I 1 *\| 

-M™ M , + C r'--g Ii L^ 



tp -log[27Vlog(t)], 

with C and to constants, so a a exp(2tp) grows slowly, 
as log(t). 

During a late-time era dominated by a constant vac- 
uum energy density, pa — 3i?p i with de Sitter expansion 
of the form a = exp(i?oi) we have late-time solutions of 
the form 

tp ~ p a + B cxp(-3g t) - N( - 3 ^ 2 + 1} cxp(-2C - 3H t), 



<P 



(12) 



where tpo, B and C are constants. This case corresponds 
to the addition of a constant potential V = Vq for the 
scalar field and we see that the effect is to turn off all time 
variations in <p, and hence a. A constant asymptote for 
tp(t) also occurs for any accelerated expansion in which 
a = t n , with n > l.Any potential of the form 

V(tp) = V Q + U(tp) 

where U falls off as exp(— ptp n ) for large tp, with n > 1, 
will result in tp approaching a constant value as t — > 00 so 
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long as the kinetic energy of the field is negligible com- 
pared to the matter density, which is the physically real- 
istic situation. In contrast, the kinetic energy \p 2 may 
dominate as t — > 0. If it does then it leads to evolution 
of the scale factor with a = i 1 / 3 and an exact solution 
for the scalar field evolution of the form [43l 



these solutions here and also derive a leading-order solu- 
tion for the scalar field in a universe dominated by dark 
energy (the dark energy is not due to the scalar field ip 
here, but to the other matter). 



<P= glogf ■jJ-log(£)+ilog(t)+log 



t to 



where E and t are constants; this solution approaches 
(p = (| - E) \og{t) as t -> and p = (± + E) log(i) as 
oo, so the fine structure constant evolves as a a 



l±2E 



in these limits if the kinetic energy dominates. 



V. THE CASE OF C = : a CONSTANT 

When £ = 0, there is no coupling of the scalar field 
to the electromagnetic matter, the problem reduces to 
the cosmology of a scalar field in the presence of a per- 
fect fluid and there is no variation of the fine struc- 
ture 'constant' a. In order to include the effects of a 
self-interaction potential, we shall assume two popular 
choices, an exponential potential and an inverse power- 
law potential for the scalar field, 



V{ip) 
V(<p) 



V exp( 

v*p-^ 



-\<p), 



in which A and 7 are dimensionless constants and Vq > 
0, 14 > are constants with dimensions [mass] 4 and 
[mass] 4+7 , respectively. These two potentials have been 
studied extensively; they can be used to obtain power-law 
inflation when the scalar field is the only matter source, 
and have scaling solutions where the energy density of 
the scalar field evolves in proportion to the density of the 
dominant fluid component of the universe in the presence 
of matter and radiation. With the potential added, the 
field equations become: 

3H 2 = Pm + p r + hp 1 + V(<p) + p&3) 

p + 3H(p + V'{p) = 0, (14) 
p m + 3Hp m = 0, (15) 

p r +AHp r = 0, (16) 



1. Radiation- dominated solution 

In the radiation-dominated era, we could neglect the 
non-relativistic matter species and vacuum energy den- 
sity (p m = pa = 0) and obtain the following solution: 



a oc t 1/2 , 
1 



H = 
Pr 



2f 



PrO 



2 t 
A l0g V 



(17) 
(18) 

(19) 
(20) 



where p r Q, to and ipo are constants. It is easy to see that 
the scalar field energy density is given by 



1 



/V = ^ + V(ip) 



(21) 



in which we have defined Vq — Vq exp(— A<^o)- Thus, 
p v cx t~ 2 scales in proportion to the radiation energy 
density p r and their ratio is kept constant during the 
evolution. 

Note the Friedmann equation and the scalar field equa- 
tion of motion give two algebraic relations between the 
constants defining the scaling solution: 



3 
4 



Pr()io + To + ^0*0 



and 1 = \ 2 Vot 2 [1 



and so 



f> 2 1 , - 2 3(A 2 -4) 
V t = — and p r0 t = — — ' 

Hence, the constant fractional energy densities are given 
by 



in which we have defined p r = p r exp(— 2ip). We shall 
discuss two particular potentials in turn below in prepa- 
ration for the discussion of the situation where varying a 
is introduced. 



A 2 -4 



A 2 



4 

A 2 ' 



(22) 
(23) 



A. Exponential potential 

For the exponential potential it is well known [3, EH 
that scaling solutions for tp exist when the universe is 
dominated by either radiation or matter. We summarize 



2. Matter-dominated solution 

Similarly, in the matter-dominated era we can neglect 
the radiation and vacuum densities to obtain the follow- 
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ing solution: 



2/3 



a (x t 



H = —, 



9 



PrnO 



2 i 
A^V 



(24) 
(25) 

(26) 
(27) 



where p m o,io and ¥>o are new constants. Here, Eq. (21) 
still holds and p v now scales as p m . The Friedmann 
equation and the scalar field equation of motion give two 
algebraic relations between these quantities 



PmO*0 + ^2 + ^0*0 



and 2 = A 2 F *o, 



the dust-dominated era. It has no closed analytical so- 
lution and we shall seek a self-consistent approximate 
solution following the logic in ref. [3l|. We start from 
the ansatz that in the A dominated era the field ip is 
slowly-rolling, and it is easy to obtain the slow-roll solu- 
tion if ~ i \og(XWt/3<;) by setting = 0. Let us next 
make the following approximation by an asymptotic se- 
ries: 



A 



log 



XW 



^ a n t~ 



n=l 



where a n are some constant coefficients. Substituting 
this back into the scalar field EOM Eq. (30) we have 

1 OO q OC 

1 i \ " / i i\ fl " i ^ \ " a n 



which lead to 

f> 2 2 2 4 (A 2 - 3) 

^0*0 = ^2 and P ,n0t = ^2 ' 

and so the constant fractional energy densities are given 
by 



A — 3 
A 2 ' 



(28) 
(29) 



Note that in order for fl v < 1 during the radiation era, 
A 2 > 4 is required, and this then ensures f2 v < 1 during 
the dust era and Cl m > 0. 



3? 
Xt 



exp 



xy ^ 



3£ 

Xt' 



as t — > oo. Choose appropriate a n so that the terms l/t r 
with r > 2 cancel, we find that the solution for ip can be 
written as 



ip - 



log 



XW 



i 



2 



V 



(31) 

It is clear that as time grows, the asymptotic series be- 
comes less important and so the slow-roll solution is ever 
improved. Eq. (31) is a good approximation when t is 
large as we have seen in the derivation; when t is small, 
Eq. (30) could be linearized as 



3. Dark- Energy -dominated solution 

Recent observations suggest that the universe is cur- 
rently, and will remain, dominated by a gravitationally 
repulsive form of matter dubbed dark energy. In the 
simplest scenario, this is just a cosmological constant for 
which the expansion rate of the Universe will tend to a 
constant, while in other models it can be exotic matter, 
or changes to the law of gravitation, which drive a differ- 
ent future evolution of the universe. Here we consider the 
evolution of our ip driven by the exponential potential in 
the background of simple dark energy domination. For 
simplicity we consider just two cases, where the cosmic 
expansion factor is either exponential, a = cxp(Hot), Hq 
constant, or a power-law in time, a oc t@ with > 1. 

a. Case 1 This is a cosmological constant (A) dom- 
inated universe, for which the scale factor evolves as 
a oc exp(^i) where <; = \/A/3 and so the EOM for ip 
becomes: 



ip + 3<rci> = W exp(— Xip) 



(30) 



where we have defined W = XVq. This equation is very 
similar to the scalar field EOM for the BSBM model in 



(p + 3^ip = W 

where we assume the initial value of ip is zero. The solu- 
tion is then 



ip = ip c + A exp(— 3?£) + 



Wt 

Ik 



Wt 
Ik 



(32) 



where (p c and A are constants of integration. The linear 
term in t seems to be the second Taylor term of the slow- 
roll solution. We can see that in both solutions the scalar 
field ip will not tend to constant when t goes large, which 
is in contrast to the case of BSBM. 

b. Case 2 This is described by a oc i" (n > 1), for 
which the scalar field EOM could be written as 



3n 

ip + — ip = W exp(— X<f) 



which has an exact solution 



A 2(Sn 



XW 



T) + -x lost 



(33) 



(34) 



Again, we find a logarithmic behaviour of if in the ac- 
celeration era, which means that ip will never approach 
a constant. This is, of course, not surprising because 
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we know that the exponential potential has tracking be- 
haviour for any power-law background expansion with 
n > 1/3, no matter whether it is n < 1 (matter and radi- 
ation dominations) or n > 1 (dark energy domination). 

Of course, we also need to justify the assumption that 
the energy density in the scalar field is always subdomi- 
nant. In case 1 this is obvious because the A has a con- 
stant energy density while the scalar held has a p/p ratio 
which is greater than —1, meaning that its energy density 
decays continuously. For case 2 we have p v cx t~ 2 cx pue 
and it again exactly tracks the dominant component in 
the universe. In both cases there is no way for tp to come 
to dominate the total energy density. 



B. Inverse power- law potential 

We turn next to the inverse power-law potential V = 
V*<p~' 1 , with 7 a positive constant. It is well known (46| 
that this potential also permits tracking behaviour of the 
scalar held tp. 

1. Radiation- and Matter- dominated solutions 

Suppose the background universe expands according to 
a oc t n and the energy density in tp is only subdominant, 
then tp has the solution: 



ip : 



(35) 



where A is constant to be fixed. To determine the value 
of A, take the time derivatives of tp 



2A 
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-t y+ 2 . 



7 + 2 
2A~f 



(7 + 2) ; 



27 + 2 
-t -T + 2 



and insert them together with H — n/t into the scalar 
held EOM, we get an algebraic equation for A which has 
the solution 



.4 = 



7(7 + 2) 2 K 
6n( 7 + 2)-2 7 



(36) 



Since tp oc 1 1+ 2 . 

2l 



it is easy to see that p v 

V(ip) oc £ - tT2 — > t~ 2 for 7 ^> 1. As the dominant com- 
ponent in the universe also has an energy density scaling 
as i~ 2 , we see that the energy density of ip simply tracks 
the dominant component, which is radiation in the radi- 
ation era and dust in matter era. 

Note that this tracking behaviour is only approxi- 
mate for 7 > 1, while in reality p v decays slower than 
^dominant- This means that the fractional energy density 
of the scalar held tp is ever increasing and eventually will 
no longer be subdominant. However, for enough large 7 
this will take a very long time so the issue will not bother 
us for some time. 



2. Dark-Energy-dominated solution 

We next consider the evolution of tp in a dark-energy- 
dominated universe. Again, we consider two cases, case 

1 for A domination where if is a constant and case 2 for 
a power-law inflation a oc t n (n > 1). Obviously case 

2 has the same behaviour as in the radiation or matter- 
dominated universes but with the value of n in Eq. (36) 
changed, and so we will not consider it again here except 
stating that in the oat™ dark energy era the held tp 
does not stop growing. 

In the first case, of A domination, we now write the 
scalar field EOM as 

7K 



tp + 3<;tp 



tpi 



The slow-roll solution to this equation is 



tp , 



7(7 + 2) V* 
3^ 



(37) 



(38) 



When t goes large, the tp term will be less and less im- 
portant because tpjtp oc 1/t and so the slow-roll solution 
is ever improved. Also, again the energy density in the 
scalar field tp decays in time so that its fractional energy 
density always decreases and it will never dominate the 
total energy density. Eq. (38) indicates that tp will con- 
tinue to grow in the A-dominated era; however the rate of 
growth is lower than that in case of a oc t™ [c.f. Eq. (35)], 
and can be very low when 7 — > 00. If needed, we could 
improve Eq. (38) by adding an asymptotic series, of which 
the leading terms are 



7+2 



7(7- 



3c 



2)K t+ 7(7 



9^ 2 



logt+const.+O 



log* 



VI. THE GENERAL CASE OF(/0 AND V ^ 

Our ultimate aim is to consider the evolution of the 
scalar held tp when both the bare potential [c.f. § 4] and 
the matter coupling term [c.f. § 3] are present, from which 
we can learn how the hne structure 'constant' a evolves in 
time. To that end we draw the contents of the above two 
sections together to get a picture of the whole evolution 
of tp in the presence of both terms. We consider again 
the two cases: an exponential potential and an inverse 
power-law potential. 



A. Exponential potential 

Before going to the general radiation and matter- 
dominated solutions for the exponential potential and a 
matter coupling term, we first consider the special case 
that arises when A = 8. The scalar field equation now 
becomes Eq. (10): 



tp + 3Htp - 



dp 



JCI 

2 — p n 

tjj 



-2tp 



(39) 
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and we still use the exponential potential given above. 

We can see that Eqs. (17 - 20) remain as in the 
radiation-dominated solution. This is easy to check be- 
cause for a oc t 1 / 2 , we have p m cx a~ 3 cx i~ 3 / 2 , e~ 2ip = 
(e~ Av ) 2/A cx i~ 4/A = t~ 1/2 , so a grows as exp(2</j) cx i 1/2 , 
and both sides of Eq. (39) scale as t~ 2 . In this special 
case the presence of the coupling term does not influence 
the overall form of the solution, although Eqs. (22-23) 
might be changed (slightly) so that the (constant) frac- 
tional energy density Q v is shifted in value. For the case 
A = 8, originally we had Q v = 1/16 when <r = 0, and in 
the radiation-dominated era p m <C p r so we expect the 
shift to be tiny. 

This discussion can be generalized to include a sub- 
sequent era dominated by a fluid with general equation 
of states (p 7^ 0). However, for a matter-dominated era, 
Eqs. (28 - 31) no longer remain an exact solution unless 
A — > oo. 

Let us now turn to the more general C 7^ cases. The 
effective total potential for the scalar field ip consists 
of two parts, the bare potential, V(ip) and the electro- 
magnetic matter couplings, which depends on a(t), via 
p m cx a~ 3 , so we can combine them as 

V eJf (<p) =V(tp) + ^ Pm exp(-2^). (40) 

The parameters used in this section are not specifically 
chosen to reproduce the observed time variation of the 
fine structure constant (which we defer to the next sec- 
tion), rather here we are concerned with the general dy- 
namics under the effective potential, which might also 
be useful for models of a scalar field with self-potential 
coupling to dark matter. 



1. Radiation- dominated era 

No matter which of the two parts to the scalar effective 
potential V e ff dominates, the field ip will grow at most 
logarithmically in the radiation era, and (except A = 8) 
the only difference between the bare-potential-dominatcd 
and the coupling-dominated solutions is the coefficient in 
front of \ogt. However, if that coefficient is very small 
then ip will remain approximately constant during the 
radiation era, as in BSBM with V = discussed above. 

As it is the radiation era, if the scalar field makes a 
significant contribution to the energy budget of the Uni- 
verse then p v ~ V{ip) ^> ^jPm exp(— 2tp), (notice also 

that — <C 1 which reduces the possible influence of the 
coupling terms even further), and thus V e ff(tp) ~ V(p)- 
On the other hand, if the scalar field constitutes only a 
very small part of the total energy density (for example, 
if A ^ 1), then V(ip) might be comparable to, or even 
much smaller than, ^j- p m exp(—2p) . 

Since the scale factor evolves as a cx t 1 / 2 whichever 
part of the effective potential dominates, we shall take 



this as the leading-order solution to the Friedmann equa- 
tion and look at the evolution of p under this condition. 
We then could rewrite the scalar field EOM as 



dx 



e^-ip(x) 
dx 



= N exp[-2cp(x)] + VFexp ( -x ) exp[-A<p(a;)](41) 



where we have defined x = logi, with N = 2^j-p m a 3 
and W = XVq constants. Clearly, the larger N is, the 
easier it is for the coupling term on the right hand side 
to dominate, and the larger W is, the easier it is for the 
potential term to dominate. 

As discussed above, we adopt a value of A > 2; the 
larger A is, the less important is the potential term for 
large ip. Because A > 2, the potential term decays faster 
than the coupling term, so if the coupling term dominates 
at initial time, the potential term will never become im- 
portant. On the other hand, if the potential term domi- 
nates at initial time, at some later time (if the radiation 
era is long enough) the effective potential will become 
dominated by the coupling term and the evolution of ip 
changes accordingly. 

These behaviours are easy to verify numerically by 
solving Eq. (41), and an example is given in Fig. [1] 



2. Matter- dominated era 

Now we turn to the matter-dominated era in which 
Pm ^S> Pr and thus the latter can be neglected. If the bare 
potential is the major part of the effective potential then, 
according to previous analysis, the scale factor scales as 
a cx t 2 / 3 ; if the coupling term dominates, Barrow et al. 
also showed HH, H3, H H El that a cx t 2 / 3 in the 
leading-order solution. So here we also assume that this 
is a good approximation in the matter-dominated era and 
look at the evolution of p on this background. 

In this case the scalar field EOM becomes 



-x d 
dx 

iVexp[- 



dx 

-2cp(x)] + VFexp (2x) exp[-X(p(x)]. (42) 



A qualitative analysis can be made as in the case of 
radiation domination, by considering the evolution with- 
out the potential or coupling term present, respectively. 
In the case where only the potential term is included, 
we have <p cx logt; if only the coupling term is presented 
then <p evolves as 1/2 log(27V logt) approximately. Now 
suppose that initially the potential term dominates over 
the coupling term, then because the former scales as t~ 2 
while the latter scales as p m exp(—2p) cx i~ 2 ~ 4 / A and 
falls faster, the potential term will become increasingly 
dominant and the coupling term will never become im- 
portant. On the other hand, if initially the coupling term 
dominates, then the potential term and the coupling term 
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FIG. 1: The evolution of ip as a function of log t in a radiation- 
dominated universe. The solid, dashed and dotted curves 
represent the total evolution, the evolution governed solely 
by the coupling term, and that governed only by the bare- 
potential term, respectively. The parameters are A = 10, 
N — 0.001 and W = 0.1; the initial conditions in the upper, 
middle and lower panels are ipi = and <pi — 0, 2. 5, and 5 
respectively. 



FIG. 2: The evolution of tp as a function of logt in a matter- 
dominated background universe. The solid, dashed and dot- 
ted curves represent the total evolution, the evolution gov- 
erned solely by the coupling term, and that governed only by 
the bare-potential term, respectively. The parameters chosen 
are A = 50, N = 1 and W = 0.1; the initial conditions in the 
upper, middle and lower panels are (pi = and ipi — 0, l,and 
2 respectively. 



scale as (logt) _A / 2 and t -2 [log(t)] -1 , respectively, and 
the former always falls off slower than the latter (how- 
ever, depending on the value of A, the dominance of the 
potential term could occur very late, much later than 
the transition to an acceleration era, so probably we will 
not see this transition during the matter epoch). Thus, 
in this case, the potential term will finally overwhelm 
the coupling term, and the full solution will then track 
the no-coupling one. These features can also be checked 



numerically. Note that increasing N or A will help the 
coupling term to dominate. The tracking is excellent in 
both regions. When the coupling term dominates the fine 
structure 'constant' evolves as a oc 2 N logt, but when the 
potential term dominates, a oc (t/to) 4 / x . 

A numerical example of the behaviours discussed above 
is shown in Fig. O 
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FIG. 3: The evolution of tp as a function of log a in a A- 
dominated universe. The solid, dashed and dotted curves rep- 
resent the total evolution, the evolution governed solely by the 
coupling term, and that governed only by the bare-potential 
term, respectively. The parameters chosen are A = 50, 
N = 0.01 and W = 0.0002; the initial conditions in the upper, 
middle and lower panels are tpi — and tfi = 0, 0.0005, and 
0.001 respectively. 



3. Dark-energy-dominated era 

A later times the universe will become dominated by 
any cosmological constant A and then the scalar field 
EOM becomes 

tp + 3^tp — Wexp(— Xtp) + Ncxp(— 3<;t) exp(—2tp). 

If the evolution is initially dominated by the bare po- 
tential term (W), then we know from the above analysis 



that tp ~ j logt. So the bare-potential term in the above 
equation scales as 1 jt while the coupling term scales ex- 
ponentially with respect to t, as t~i exp(— 3<;t). For large 
t, the latter decays faster, and the bare-potential term 
will eventually dominate over the coupling term and the 
fine structure 'constant' will evolve all the time in the fu- 
ture. A numerical example is shown in Fig. [3l Note that 
for this analysis we have defined a different set of param- 
eters x = logo, with W = AVo/c 2 and N = 2\(\p m o/<; 2 
constants, so that the above EOM is rewritten as 



•3^-¥>(£c) 
ax 



= N exp(-3x) exp(-2tp) + Wexp(-Xtp). 

If the dark energy drives power-law inflation a cx t n 
(n > 1) of the universe today, then the analysis of 
the whole evolution is qualitatively the same as for a 
radiation-dominated universe and will not be repeated 
here. 



B. Inverse power- law potential 

Next, we consider the evolution of tp under the con- 
trols of the coupling term and a bare inverse power-law 
potential. 



1. Radiation- dominated era 

In the radiation-dominated era the cosmic scale factor 
scales as a oc t 1 / 2 , while the energy densities of scalar 
field and dust can be neglected. So, just as in the case of 
exponential potential, we can write the scalar field EOM 
as 



,d_ 

dx 



e ~' x Tx^ x) 



= N exp[-2ip(x)] + Wexp i-xj <p~ (l+1 \ (43) 

where again x = log t, but now W = jV* . 

According the above analysis, if the bare poten- 

2 

tial term dominates then <p evolves as ip ~ t~/+ 2 ~ 
exp ^^5^^ j (so tp is exponential in x), while if the cou- 
pling term dominates then ip ~ j^ogt ~ ix (so <p is 
linear in x). These features can be seen clearly in the nu- 
merical example given in Fig. 31 Note that in the scalar 
field EOM the bare-potential term decays as ~ <p~^ +1 \ 
while the coupling term scales either as ~ p m exp(— 2ip) ~ 
t~z exp(— 2ip) ~ ( y 9~4(T+ 2 ) exp(— 2tp), (in case the effec- 
tive potential V e f / is dominated by the bare potential) , or 
as ~ p m exp(— 2ip) ~ exp(— 2ip) ~ exp(— 8^), (when 
V e ff is dominated by the coupling term), so obviously in 
both cases when eventually tp is large enough the bare- 
potential term will dominate over the coupling term driv- 
ing the evolution of tp. This can also be seen in the figure. 
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FIG. 4: The evolution of ip as a function of log t in a radiation- 
dominated universe. The solid, dashed, and dotted curves 
represent the total evolution, the evolution governed solely 
by the coupling term, and that governed only by the bare- 
potential term, respectively. The parameters used are 7 = 50, 
N = 0.001 and W = 0.005; the initial conditions in the upper, 
middle and lower panels are (pi = and tpi = 1,5, and 10 
respectively. 



FIG. 5: The evolution of tp as a function of logt in a matter- 
dominated universe. The solid, dashed and dotted curves 
represent the total evolution, the evolution governed solely 
by the coupling term, and that governed only by the bare- 
potential term, respectively. The parameters chosen are A = 
100, N = 10 and W = 0.005; the initial conditions in the 
upper, middle and lower panels are tpi — and <pt = 1, 2, and 
3 respectively. 



EOM becomes 



2. Matter-dominated era 



In the matter-dominated era the scale factor evolves 
as a oc i 2 / 3 , while the energy densities of scalar field and 
radiation can be neglected. In this case the scalar field 



dx 



= N exp[-2<p(x)] + Wexp(2x)(p~ h+1 \ (44) 

where x and W are as defined above. 

From the above analysis we know that if the bare po- 

2 

tential term dominates then ip evolves as <p ~ t~>+ 2 ~ 
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FIG. 6: The evolution of tp as a function of log a in a A- 
dominated universe. The solid, dashed, and dotted curves 
represent the total evolution, the evolution governed solely 
by the coupling term, and that governed only by the bare- 
potential term respectively. The parameters are 7 = 50, N = 
W = 0.01, and the initial conditions are tpi — 0, and tpi = 1.1. 



exp {^pp;^ (i.e., tp is exponential in x), while if the 
coupling term dominates then to the leading order (p ~ 
ilog(logi) ~ ^logx (i.e., tp is logarithmic in x). We 
can also see these behaviours clearly in the numerical re- 
sults plotted in Fig. In the scalar field EOM the bare- 
potential term decays as ~ p~^ +1 \ while the coupling 
term scales either as ~ p m exp(—2tp) ~ t~ 2 exp(— 2ip) 



-(7+2) 



cxp(— 2tp), (when V e ff is dominated by the bare 



potential), or as ~ p m exp(— 2tp) ~ t~ 2 exp(— 2ip) 
exp [— 2 exp(2(y9)] exp(— 2tp) (when V e ff is dominated by 
the coupling). Therefore, when tp eventually grows large 
enough, the bare-potential term will dominate over the 
coupling term and drive the evolution of tp. This is ver- 
ified in the numerical results, where we can see specific 
tracking solutions in different epochs. 



3. Dark-energy-dominated era 

In a universe dominated by a cosmological constant, 
the field equation for tp in the case of a power law poten- 
tial can be written as 



d 2 d . . 

—tp{x) +3—tp{x) = 
dx z dx 



W 



tpi 



j + N exp(— 3a;) exp(— 2tp), 



where x = loga, W = 7V,/? 2 and N = 2\(\p m0 /<; 2 - 

According to the above analysis, when the effective 
potential V e f f is dominated by the coupling term, tp will 
approach a constant in the A-dominated era, while if V e f / 
is dominated by the bare potential then tp evolves as tp ~ 



3.000003 




log(a) 




log(a) 



FIG. 7; Comparison of the entire cosmological evolution of tp 
as a function of log a in BSBM model (dashed curve) and the 
model with an exponential potential V = Vbexp( — Xtp) plus 
coupling term (solid curve). The parameters for the upper 
panel are, respectively, A = 5 x 10 6 , |C|/ W — 10~ 4 and A = 
5 x 10 6 , \C\/t° = 2 x 10~ 4 . In both cases the initial conditions 
are tp i = 3,tpi = 0. 



0: (log a) ^+2. These qualitative behaviours can be 
observed in Fig. As tp goes large, the bare potential 
term in V e ff decreases as ~ tp^^ +1 '> while the coupling 
term decays as ~ exp(— 3?t) exp(— 2tp), so eventually the 
former will always dominate over the latter, driving the 
continuous growth of tp in contrast to the asymptotically 
constant behaviour seen in BSBM where W = 0. 



C. Summary of the behaviours of tp(t) and a(t) 

We summarise the possible behaviours of tp(t) found 
in the different varying-a models that we have discussed, 
in the three different cosmic eras and different situa- 
tions (bare potential V e ff -dominated or coupling £- 
dominated) in the Table shown here. The time evolution 
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TABLE I: Tracking behaviours of a in different varying-Q models in different limits. V denotes the bare potential term and C 
is the coupling term; — > means asymptotic approach. 



Model 



Radiation Era 



Matter Era 



— r - 

1 2 



Acceleration Era 



BSBM 
BSBM + Voexp(-A^) 

BSBM + V»V" 7 



exp 
i 



t * , V dom. 
£2 , C dom. 

(2t^\ , V dom. 

C dom. 



q ~ log t 
ti , V dom. 
logt, C dom. 

exp (2t^^j , V dom. 
logt, C dom. 



a - 

4 

«*, 

2 



const, 
a oc t n ; 
H = const. 



exp 



(2t^h) , aat"; 
exp (^2£^T2^ , # = const. 



of a is obtained from that of if by using a oc exp(2y>). 
Recall from the figures above that the transitions from 
bare-potential domination to coupling domination (or 
vice versa) are smooth, so there is a simple pattern for 
the overall ip evolution. However, note that depending on 
the initial conditions for ip, the above 'tracking' solutions 
may only be reached after a long time (as can be seen in 
the figures). In this analysis we have focussed upon ex- 
tracting the evolution of <p(t), and hence of a(t), but a 
more detailed study with a different emphasis could also 
seek the best-fit observational parameter set when vary- 
ing alpha is included, as was done in Ref. [47| . 



VII. NUMERICAL EXAMPLES OF THE ip AND 
a EVOLUTION 

We shall now consider the numerical evolution of if. 
and a, through the entire cosmological history, and try 
to connect this evolution to the observations constrain- 
ing possible time variation in a. As we have seen above, 
the evolution in if is controlled by the competition be- 
tween the coupling term and bare-potential term in the 
effective potential. The parameter |C|/u> determines the 
strength of coupling and so increasing it will increase the 
rate of variation of if. In addition, in the case of the 
exponential potential, V = V^e exp(— \<p), the parame- 
ter A controls the scaling solution of <p, and the larger 
it is, the smaller the fraction of the total energy density 
Sip tends to be; note also that A governs the slope of 
the evolution of if : if there is no coupling term then 
the solution of a oc exp(2ip) is given by a oc a - a in the 
radiation era and a oc a~~ in the radiation era. This 
power-law evolution means that in order that the frac- 
tional variation of a between now and z ~ 6 should not 
exceed the observational bounds, i.e., Aa/a < 0(1O -5 ), 
A must be very large. Similarly, in the case of an inverse 
power-law potential, 7 controls the scaling solution of if, 
and larger values of 7 correspond to smaller variations of 
ip in time, so that an observational constraint that the 
allowed variation of a oc exp(2(p) be small requires 7 to 
be very large. 

Finally, the initial value ifi, is also an important quan- 
tity; although for the exponential potential one can al- 
ways choose ifi = by adjusting Vq correspondingly. For 



the coupling term we do not have this freedom - a larger 
ifi will weaken the coupling through exp(— 2<fi) and thus 
reduce the change in ip. 

In Figs. [7] and [5] we have plotted the evolution of f{t) 
for some choices of the model parameters for the two 
above-mentioned potentials, respectively. For compar- 
ison, we also plot the result for BSBM model (where 
V = 0) with the same choice of the parameter |C|/<x>. 
In both cases <p is initially dominated by the bare- 
potential term but later becomes dominated by the cou- 
pling term because the largeness of A (or 7) means that 
the slope of the bare-potential-dominatcd evolution is 
smaller that the coupling-dominated evolution at later 
matter-dominated times. When the universe becomes 
dominated by the cosmological constant, the coupling- 
dominated solution for ip. approaches a constant and the 
bare-potential-dominated solution grows very slowly, so 
the total solution grows very slowly too (almost con- 
stant in the figure). Note that the addition of a bare- 
potential term makes if begin evolving earlier. This pro- 
duces an earlier onset of variation for a than in the pure 
BSBM model (the late-time evolution, which is relevant 
for the QSO observations, is however almost the same as 
in BSBM because at this late stage the total solution is 
dominated by the coupling term). 

Figures [9] and \W\ show the evolution of a/ao (where 
ctQ is the current value of the fine structure constant) in 
the two models when compared with the prediction of 
the BSBM model (dashed curves). The qualitative fea- 
ture in a model with self-potential for if is the same as 
in BSBM: at early times a remains a constant; during 
the matter-dominated era there is a slow growth; and 
then, when the cosmological constant begins to domi- 
nate, the growth stops. The differences are: first, the 
commencement of growth for a starts earlier than the 
BSBM case because in this model the bare-potential term 
could drive the evolution of if even in the radiation- 
dominated era; second, the late-time evolution is dom- 
inated by the coupling term and so the late-time evo- 
lution of if is similar to that in BSBM, but the ear- 
lier commencement of the growth means that the to- 
tal variation of a is greater than that in BSBM. The 
fractional change of a between z ~ 4 and now is about 
0.5 x 10 -5 which is consistent with the QSO observa- 
tions. The current rate of tp is given by ipo = ( •£ ) Hq 
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FIG. 8: Comparison of the entire cosmological evolution of ip 
as a function of log a in the BSBM model (dashed curve) 
and the model with an inverse power-law potential V = 
V*ip~~" plus coupling term (solid curve). The parameters for 
the upper panel are, respectively, 7 = 5x 10 6 , \(\/lo = 2x 10~ 6 
and 7 = 1 x 10 7 , \(\/lu = 2 x 1(T 6 . In both cases the initial 
conditions are tpi = 1, tpi = 0. 
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FIG. 9: The evolution of a/an versus log a in the BSBM 
model (dashed curve) and in the model with an exponential 
potential plus coupling term (solid curve). The parameters 
chosen are A = 5 x 10 6 , [C|/o> = 10~ 4 . 




where Hq w 70 kms 1 Mpc 1 is the present Hubble ex- 
pansion rate. For the parameters in Figure 12, we have 

(if) ~ °- 8 x 10_6 and so we have <Po ~ 0.6 x 10" 16 /yr 
which leads to (a/a)o ~ 2ip ~ 1.2 x 10 _16 /yr. This rate 
is well within all old limits Q , 0] , [H , @ but is about an 
order of magnitude above the proposed new upper bound 
(l9| on the current rate of a variation from atomic clocks, 
which is (d/a) = (-1.6 ±2.3) x lCr 17 /yr, although the 
uncertainties may be modulated slightly by accounting 
for seasonal variations in the local gravitational poten- 
tial 0. 



FIG. 10: The evolution of a/an versus log a in the BSBM 
model (dashed curve) and in the model with an inverse power- 
law potential plus coupling term (solid curve). The parame- 
ters chosen are 7 = 5 x 10 6 , [C[/a> = 2x 10~ 6 . 



VIII. SUMMARY AND CONCLUSIONS 

In this paper we have considered the dynamics of the 
varying-a theories which arise when the original BSBM 
theory is generalised by introducing an exponential or 
inverse power-law self-potential for the scalar field driv- 
ing the variation of a. These two representative poten- 
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tials capture the essential ingredients of general poten- 
tials without minima. There are two situations to dis- 
tinguish and analyse separately: according as to whether 
or not the scalar-field potential comes to dominate the 
late-time dynamics of the universe. In combination with 
the coupling with matter, the additional bare potential 
forms an effective total potential V e f / for the scalar field 
ip which governs the allowed time variation of a. We 
have presented the solutions to the scalar-field equation 
of motion in cases where V e fj is dominated solely by 
the coupling term or the bare potential, respectively, in 
different cosmic eras. In most cases the bare-potential- 
dominated solution differs from the coupling-dominated 
one; the contributions of these two terms to V e ff vary 
with time, and it is possible for there to be a transition 
from one solution-type to another. The numerical results 
show that the transition between solutions types can be 
very smooth, and for most of the time the solution for the 
scalar field tracks either the bare-potential-dominated or 
the coupling-dominated solution. These features ensure 
that the evolution of <p under V e f f has a very simple pat- 
tern. The main results are briefly summarized in Table 
1, and these results are quite general, not depending on 
whether the parameters defining the potential (A and 7) 
are extremely large or not. 

The consequences for the time evolution of the fine 
structure 'constant' of adding potentials V(tp) to the 
BSBM theory are summarised as follows. In light of 
the observational constraints on how much variation in a 
there can be over redshifts z < 6, we find that the restric- 



tions on the interaction potential parameters (A in the ex- 
ponential potential and 7 in the inverse power-law poten- 
tial) must be very strong, in order to prevent the bare po- 
tentials from becoming unacceptably dominant and driv- 
ing unacceptable fast time variation of a. For example, in 
the case of an exponential potential, V — Vq exp(— Aip), 
we must have A > 10 6 - 10 7 with \C\/w ~ 0(l(r 6 ), and 
for the inverse power-law potential, V = V*(p~ 7 , we need 
7 > 10 6 - 10 7 with \C\/oj ~ 0(l(r 6 ) - the exact con- 
straint, of course, depends also on the initial conditions). 
This means that the A and 7 are constrained to be so 
large that if they appear in quintessence models then the 
scalar field is practically indistinguishable from a cosmo- 
logical constant, which has no dynamics at all. The total 
variation of a can be enhanced compared to the case of 
no bare potential (BSBM), and the variation commences 
much earlier. Finally, because with an exponential or 
inverse power-law potential the scalar field will not ap- 
proach a constant even in the dark-energy-dominated era, 
the fine structure constant a will continue to increase for 
all future time, until eventually there will no stable atoms 
in the Universe at all [32| . 
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